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Abstract 

We describe methods for evaluating one-loop integrals in 4— 2e dimensions. We give a recursion 
relation that expresses the scalar n-point integral as a cyclicly symmetric combination of (n — 1)- 
point integrals. The computation of such integrals thus reduces to the calculation of box diagrams 
(n = 4) . The tensor integrals required in gauge theory may be obtained by differentiating the scalar 
integral with respect to certain combinations of the kinematic variables. Such relations also lead 
to differential equations for scalar integrals. For box integrals with massless internal lines these 
differential equations are easy to solve. 
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Many processes of interest at current and future e + e~ and hadron colliders involve large 
numbers of final state particles. Radiative corrections to these processes are needed for precise tests 
of the standard model. It is therefore useful to have techniques for evaluating one-loop integrals 
where the number of external legs is large. As the loop integrals appearing in radiative corrections 
are often infrared and/or ultraviolet divergent, it is desirable to regulate them by performing them 
in 4 — 2e dimensions. 

In this Letter, we derive a relation between the n-point and (n — l)-point one-loop integrals, 
which for n > 4 allows the recursive determination of the general n-point scalar integral in D = 
4 — 2e, as a linear combination of box integrals (n = 4), provided only that the external momenta are 
restricted to lie in four dimensions, and neglecting 0{e) corrections. The required box integrals can 
generally be evaluated in closed form through 0(1), that is to say, with all poles in e manifest, and 
with all functions of the kinematic invariants expressed in terms of logarithms and dilogarithms. (A 
compact expression for the general infrared-finite box integral has recently been given by Denner, 
Nierste, and Scharf [1]; the infrared-divergent box integrals with all internal lines massless are 
collected in ref. [2].) Therefore, the higher-point integrals can now be represented in the same 
closed form. In a separate paper [2], we apply these techniques to determine explicitly the pentagon 
integral with all external lines massless, or with one external mass. 

Various authors [3,4,5,6] have discussed the computation of pentagon and higher-point integrals 
that can be evaluated in D = 4 (i.e. that are infrared finite). In particular, Melrose [3] and 
independently van Neerven and Vermaseren [5] have expressed the D = 4 pentagon integral as a 
linear combination of five D = 4 box integrals, and the relation we find for n = 5 may be thought 
of as the dimensionally-regulated version of their equations. References [3,5] also express the four- 
dimensional n-point scalar integral for n > 6 (with external momenta restricted to D = 4) as a 
sum of six (n — l)-point integrals; the derivation in ref. [5] extends straightforwardly to (4 — 2e)- 
dimensional loop-momenta as well. For n > 6 these relations are of a somewhat different type than 
the relations that we find. We have been informed that Ellis, Giele, and Yehudai [7] have recently 
evaluated the D = 4 — 2e pentagon integrals by an independent technique. 

In gauge theories, tensor integrals appear in which the n-point integral may contain up to 
n powers of the loop momentum in the numerator of the integrand. It is possible to perform a 
Brown-Feynman [8] or Passarino-Veltman [9] reduction of the integrand, solving a system of al- 
gebraic equations to reduce the tensor integrals to a linear combination of scalar integrals [10]. 
The framework developed here provides another method for computing tensor integrals. Feynman 
parametrization converts tensor integrals into integrals where polynomials in the Feynman param- 
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eters (of order < n) are inserted into the numerator of the integrand. (In the string-motivated 
technique [11] for evaluating QCD amplitudes such a representation is obtained directly.) We 
will see that such integrals may be obtained by differentiating the scalar integral with respect to 
particular combinations of the kinematic variables. 

The basic integral we wish to evaluate is the dimensionally-regulated one-loop scalar integral 
with n external momenta ki, n external masses rrtj (kf = mf), and n internal masses Mj, i = 
1,2,. ..,n: 

2e f 



In = A* 

where we take 



(2vr )4- 2e (p 2 _ M 2) ( (p _ pi) 2 _ M 2) . . . ((p _ pn _ l? _ M 2) > 



PO = Pn = 0. 



(1) 



(2) 



Introducing Feynman parameters in (1), J n becomes 

■1 , d 4-2e p 



l n = n e (n - 



n J_ n 7^. ^ £ -IV 



(3) 



Completing the square in the denominator, Wick rotating, and integrating out the loop momentum 
yields 



I n [l] = T(n-2 + e 



[ tf 
Jo 



di 5(1 - Y,i a i) 



n-2+e 



where the scalar denominator V(a.i) is 

V{ ai ) 

and we have defined 



L i=l 



1 2 



i=l 



/„[!] = i (-1)" +1 (47r) 2_e ^ _2e X„ 



(4) 



(5) 



(6) 



as a convenient normalization for the Feynman parametrized scalar integral. (The expression inside 
the brackets indicates the Feynman-parameter polynomial in the numerator of the integrand; non- 
trivial polynomials correspond to tensor integrals.) We rewrite the second set of terms in P(oj) 
using J2i a» = 1, in order to make it homogeneous of degree two in the a, and symmetric in a, «-> af 



where the matrix Sij is given by 



23(oi) = S i 



Stj = \{M? + M*-p%). 
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(7) 



(8) 



Here pij is the sum of \i — j\ adjacent momenta, 

Pa = 0, = pj- X -Vi-\ = ki + k i+1 -\ Vkj-i iox i < j. (9) 

Now we are in a position to derive some useful general relations. We start by evaluating the 
following integral two different ways: 

j = , r _ n n+l U ^2-e f d*- 2 * P p 2 ~Ml 

11 j [ ' J (2vr) 4 "^ (p2 - Ml) {{p - Pl y - Ml) • • • ((p - p n _i) 2 - M 2 ) " UUj 
First, we cancel the numerator against the denominator, yielding the (n — l)-point integral 

Jl = 'till [1] • (11) 

In general, the notation i^A^l] refers to an (n — l)-point integral whose kinematics is 'inherited' 
from its 'parent' i n [l] by removing the propagator between legs i — 1 and i. Second, we evaluate 
Ji as an n-point integral. Feynman parametrizing the integral and completing the square in the 
denominator in the usual way yields 

d^q 

1 ' (12) 



Ji = i(-l) n+1 (4vr) 2 - e (n-l)! J * d n a, 5(1 - | 

9 2 +^+Er=i^(pLi-M 2 -M 2 ) 



[g 2 - P]« 

Next we integrate out the loop momentum in eq. (12) (after a Wick rotation), and equate the result 
to the (n — l)-point result (11), to get an equation for I n [a,i]: 

n 

-I { n\[l] = [-i(4-2e) + (n-3 + e)]/f= 6 - 2e [l] + £>?_!- - M?)J»[ai] 

(13) 

= ( n -5 + 2e)/,f= 6 - 2e [l] - 2^S H /„N • 



i=l 

,2 



In this equation we have rewritten the terms coming from q +T> in the numerator of equation (12), 
using the Feynman-parameter representation of the scalar n-point integral in D = 6 — 2e dimensions, 



/f= 6 - 2£ [l] = r(n - 3 + e) jf 1 d% <5(1 - £,«,) 



n— 3+e 



(14) 



(This representation can be obtained from equation (4) simply by letting e — ► e — 1, which shifts 
D = A-2etoD = 6- 2e.) The properties of 7 J f =6 - 2e [l] as e will play a role below. 

Similarly, by considering the integral J, with (p — Pi-i) 2 — M 2 in the numerator, we find the 
set of equations 

n 

2 /„[«,-] = &[!] + (n-5 + 2 e )/,f= 6 - 2 ni] . (15) 



Solving these equations for I n [ai] we get 



In M 



where 



E^tf-iM + (n - 5 + 2e) * 



E s «'- 

i=i 



= 6-2e 



[1] 



(16) 



(17) 



Now sum equation (16) over i and use ^E a» = 1 to get 



InN 



.i=l 



Ci + (n-5 + 2e)c I 



D=6-2e 



[1] 



where 



Co 



i=i 

r D=6-2e 



E V- 



(18) 



(19) 



Finally, use equation (18) to eliminate J 5 6 from equation (16), and thereby obtain an equation 
for the one-parameter integrals I n [a,i] in terms of I n _i[l] and /„[!]: 



In[di 



1 ™ 



where 



j=i 



. | Ci£j 

'•J — a ij 



CO 



(20) 



(21) 



The external momenta as well as the loop-momenta in equations (16), (18) and (20) are still 
(4 — 2e)-dimensional. In computing one-loop corrections to physical processes, one is generally 
interested in restricting the external momenta to D = 4. On the other hand, for n > 6 the 
coefficients Cq and q given in equations (19) and (17) appear to be singular for D = 4 kinematics. 
(The rank of the n x n matrix Sij is n — 6 in D = 4 [3], so for n > 6, S is not invertible.) In fact cq 
and Ci are nonsingular in the limit of D = 4 kinematics (co actually vanishes for n > 5). In order 
to see this, and in order to see how to apply these results to tensor integrals, it is useful to perform 
two changes of variables: first a change of integration variables in the integral (4), then a change 
of kinematic variables. 

Following 't Hooft and Veltman [4] , we make the change of integration variables in equation (4) , 



a* = 



a n = 



CtiUi 



En 
j=l a 3 U 3 



no sum on i, 



En— 1 
3 = 1 U 



(22) 



En 
3 = 1 a 3 



U-j 



Assuming that all ccj are real and positive (physical regions may be obtained by analytic continu- 
ation), the integral becomes 

11^ / d ^ f — T^+l . (23) 

j=i / [E" J= i/w% 

where 

PlJ = Sij-Qiaj = \(Mf + M 2 -p 2 j )a l a ] . (24) 

The form (23) for the integral will be the most useful if we also make a change of kinematic 
variables such that the scalar denominator (i.e. the matrix p) no longer depends on the n variables 
«j. First we describe a linearly independent set of the original kinematic variables. The n(n + 
l)/2 Mandelstam variables = (A;,- + kj) 2 are not linearly independent, due to n momentum 
conservation relations. Instead one may use the square of the sum of I adjacent momenta, p 2 i+l = 
{ki + ki+i + . . . + ki + £-i) 2 . For n odd, t runs from 1 to (n — l)/2. For n even, t runs from 1 to 
n/2, but for I = n/2 there are only n/2, rather than n, such variables. These n(n — l)/2 linear 
combinations of Mandelstam variables form a linearly independent, cyclicly symmetric set. We 
should add to this set the n internal masses M 2 . (Alternatively one may use the coefficients Si^ + e 
and Su appearing in the scalar denominator P(aj).) For example, for the hexagon (n = 6) the 
independent variables are M 2 , mf, Sj,j+i = (ki + fci+i) 2 , for i = 1, 2, 3, 4, 5, 6 (all indices are taken 
mod 6), along with t M+M+2 = (A;,- + k i+1 + /c i+2 ) 2 , for i = 1, 2, 3. 

The following is an example of a change of kinematic variables that eliminates all a^-dependence 
from the scalar denominator: 

{M;, 7713,^+1,^+1,3+2, • • •} -> {ai]M 2 ,m 2 ,ii, . . .} , (25) 

where the new variables {c^; M 2 ,m 2 ,ii, . . .} are defined by 



T77 2 









777 2 


M 2 




ttiOi+1 


~ «? 


«?+i 


1 


M 2 




"iOi+2 


" «? 


a i+2 




M 2 


Mf+3 


aiOi+3 


«? 


«?+3 



(26) 



*i,i+l,i+2 — 

«i«i+3 «j "i+3 

(i = 1,2, ... ,n), 

with all indices taken mod n. This is certainly not the only change of variables possible, and in 
some cases it is convenient to make other choices. (Indeed, for the special cases n = 8, 12, 16, . . ., 
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equations (26) are not a legitimate change of variables. This problem can be cured by modifying 
slightly the substitution for just two of the s^+i, say S12 = ~ ~ an d s 23 



Mi 



Mi 



so that the set of kinematic variables is now {af, Ai, A2, Mf,rhf,ti, ■■■}■) 



Q2Q4 C*2 

With this change of kinematic variables, tensor integrals may be calculated simply by differ- 
entiating the scalar integral with respect to the a^. Consider the Feynman-parametrized integral 
with an arbitrary monomial of degree m inserted in the numerator, 



Inia^a^ ...a im ] = T(n-2 + e) / d n ai 5(1 - Y.i a i) 



s: 

Jo 



(Xj, Q>i 



a; 



[V{oi 



ln-2+e 



(27) 



Using equation (23) it can be represented as 

T(n-3-m + 2e) ( A 
I n [a tl a l2 ...a lm ] - — _ - - — m, 

If we define the reduced integrals 



_d d f I n [l] 



• (28) 



L = 



Haj) I^Pdai/ai})] , 
7=1 



= (n«i) . 



(29) 



then we may write 



In \o>i\ ®i2 • • • a i 



r(n - 3 - m + 2e) d m /„ 



(30) 



T(n-3 + 2e) da h ■ ■ ■ 8a im ' 

We now proceed to find simple expressions for the coefficients cq, Q and Cjj appearing in 
equations (18), (20), and (16), in terms of the a, variables. To do this it is first useful to examine 
how the Gram determinant of the (n — l)-vector system associated with the n-point process depends 
on «j. The Gram determinant is defined by 



A n = det (2ki ■ kj), 



(31) 



where the prime signifies that one of the n vectors ki is to be omitted before taking the determinant; 
due to momentum conservation, ^ ki = 0, any one of the vectors may be omitted. ^ 
Next we introduce the reseated Gram determinant, 



n 

A n = (n«*) A? 



(32) 



t The notation for, and normalization of, the Gram determinant in equation (31) differ from other conventions 
in the literature, e.g. references [3,12]. 



which will turn out to have several useful properties. Rewrite A n in terms of the variables and 
the matrix p defined in equation (24), after omitting k n in the definition of A n : 

n n 

An = (T\aj) det (2ki ■ kj) = (t[ a i) de } " Pi) 

e=i e=i 

= 2 n - l (f[ai) detf-^---^---^- + ^ (33) 

\- LJ - / \OiiOij OLiOL\ OL\OLj a.\ ) 



.1... f .. .. a 3 .. a i a 3 



= 2 n det \ Pij - Pi\— - p\j h Pn 

«i «i «i 

By omitting other vectors fcj in the definition of A n , it is easy to obtain n different expressions for 
the same quantity A n , 

A„ = 2"" 1 a, 2 det -p*A + ) , £=l,2,...,n. (34) 

i,3^t\ on at a\ J 

From equation (34) , and the fact that the matrix p is independent of the on , it is clear that A n is 
homogeneous of degree 2 in the ctj, and that no a, appears with a negative power in A n . So we 
may write 

n 

A n = ^ Vij^iaj, (35) 

where t?^ is independent of the ctj. 

We can relate the two matrices rj and p. Consider first the diagonal element rju- The a\ term 
in A n comes from taking only the pij terms in (34), and is given by 

r] U aj = 2 n ~ l det ( Pij ) aj = 2"" 1 det p (p~ 1 ) ll aj (36) 

(where the second determinant is over all indices). Similarly, we pick off the aga m terms (m ^ £) 
in (34), by using a p^ or p ej term in place of p im or p mj -, to get r] £m = 2 n ~ 1 det p (p -1 )^. Thus 
r/ is proportional to the inverse of p, 

P = iVn rT\ r? = N n p~\ (37) 

where the proportionality constant is N n = 2 n_1 det p. 

Now we rewrite the coefficients cq, Ci and in terms of the cc, variables. First of all, using 
equations (24) and (37) the matrix S" -1 is given by 

= «< Pij «i = — ^ • (38) 

Summing this equation over i and/or j, and using equation (35) for A n , plus the definition 

ISA 



7i 



2 9a 



(39) 

non-ai variables fixed 



we find that 



cn = — c = c- = j f r? £122. ^ C40) 

In terms of the reduced integrals defined in equation (29) — which have simple differentiation 
properties — equations (18), (20), and (16) also take on a simpler form. It is convenient to make 
the change of kinematic variables for the integrals /^,-Jl], such that the variables on are identical 
to those for the parent integral I n [l]. Then we have 



Ir, = 



1 

2Nl 



X>i£li + (n-5 + 2e)A n /, 



D=6-2e 



i=l 



1 



n 



1 



4 + 2e 5a 



<9/„ f 



1 

2K 



dl„ 



n 



4 + 2e <9a; 



In[ai 



1 

2A^ 



i=i 

fU) 



n-l 



+ 



tD=6-2e 



(41) 
(42) 
(43) 



J>i#-i + (n-5 + 2e) 7i 
-i=i 

The simple equations (41), (42), and (43) are in some sense the main results of this letter. All 
the equations derived so far are valid in an arbitrary spacetime dimension D = 4 — 2e; that is, we 
have not yet assumed that e is small, nor that the external momenta lie in D = 4. However, the 
equations have the most utility in the context of dimensional regularization, i.e. for D = 4 — 2e 
with e tending to zero. We now discuss how to use the equations to recursively generate one-loop 
n-point integrals, up to 0(e) corrections. 

At first sight equation (41) does not appear very useful, due to the presence of the integral 
I^ =e ~ 2e [l] on the right-hand-side. However, for n = 5 the coefficient of this term is of order 
e; and the integral I®~ 6 ~ 2e [l] is finite as e — ► 0, because the D = 6 scalar pentagon integral 
possesses neither ultraviolet divergences nor infrared divergences (soft or collinear). So to order e 
the D = 4 — 2e scalar pentagon is given by a sum of five scalar boxes, 



5 



a l7i 4 l) [l] + 0(e). 



2N 5 



(44) 



=i 



It is easy to check that the coefficients of the box integrals in this equation are identical to those 
in the corresponding D = 4 relation in ref. [3]; both sets of coefficients are expressed in terms of 
matrix elements of Sij = \(M? + M? — pfj). In ref. [5] the box coefficients in the D = 4 relation 
were expressed using Levi-Civita symbols; in ref. [2] we show that the coefficients are nevertheless 
the same as those in equation (44). 

For n > 6, the coefficient of I® =6 ~ 2e [l] in equation (18) is not of order e. On the other hand, 
if all external momenta are chosen to lie in D = 4 dimensions, then we can set Co = 0. This is 



because the Gram determinant A n vanishes in D = 4, due to the linear dependence of the (n — 1) 
vectors k\, /C2, . . . , fc n _i [3,12,13]. In next-to-leading-order calculations, it is always possible to put 
this restriction on the external kinematics. We then get 



For n = 6 when all integrals appearing in (45) are finite, this equation is identical to the non 
dimensionally-regulated result of ref. [3], and it can be shown to be equivalent to the corresponding 
result of ref. [5] , expressed in terms of Levi-Civita symbols; the latter derivation extends straight- 
forwardly to (4 — 2e)-dimensional loop momenta. 

For n > 6, equation (45) differs somewhat from the results of references [3,5] in that it contains 
n, rather than six, (n — l)-point integrals. In the reductions (44) and (45), the cyclic symmetry of 
I n [l] is kept manifest. A more important difference arises if one wishes to extract tensor integrals 
via the differentiation formula (30). Namely, one should not restrict to four-dimensional external 
kinematics until after carrying out the differentiation, and so the cyclicly symmetric representa- 
tion (41), with unrestricted kinematics, should be used as the starting point. 

As mentioned above, for n > 6 the representations (19) and (17) of the coefficients cq and Cj 
in terms of the original kinematic variables Sij are problematic when the external momenta are 
restricted to D = 4. This is because the rank of the n x n matrix is n — 6 in D = 4 [3], so for 
n > 6 the inverse does not exist. On the other hand, the o^-representations of Co and q in 
equation (40) are non-singular and well-defined for D = 4 kinematics. In summary, the combination 
of equations (44) and (45) recursively determines the general one-loop n-point scalar integral in 
D = 4 — 2e dimensions as a linear combination of box integrals. 

Equation (18) also has significance for n < 4, even though the term containing I,f =6_2e [l] in the 
equation may no longer be neglected. For example, for n = 4 the decomposition (18) of if =4_2e [l] 
has the virtue of putting all the e — > divergences into the triangle integrals 1^ [1], since the D = 6 
scalar box is infrared and ultraviolet finite. This is important for practical calculations because 
infrared-divergent triangle integrals are generally simple to evaluate analytically; the infrared-finite 
triangles and the D = 6 scalar box can be evaluated numerically if necessary. 

Equation (42) gives a set of partial differential equations for n-point scalar integrals. For box 
integrals (n = 4), one may also use equation (43), plus the finiteness of the D = 6 scalar box, to 
simplify the differential equations through 0(e): 




n > 6, D = 4 external kinematics. 



(45) 



doLi 



4 



n 



+ (-l + 2e) 7i It 



D=6-2e 



6 



+ O(e) . 



(46) 
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In ref. [2] we solve these differential equations for box integrals with all internal lines massless, but 
with nonzero masses for 0, 1, 2, or 3 external lines. (Some of these integrals have been computed 
previously by other techniques.) Together with the infrared- finite box integral with all four external 
lines massive [1], these are the complete set of box integrals needed to recursively determine the 
n-point integrals for next-to-leading-order calculations in QCD with massless quarks. 

For example, the box integral with one external mass is needed to obtain the all-massless 
pentagon via equation (44). The solution to the differential equations for the reduced box integral 
is [2] 



/f = 2r r 



(a - +2 "- 2)£ +Li 2 (l-^)+U 2 (l-^V- 



+ 0(e), (47) 



where tt = T(l + e)T 2 (l — e)/T(l — 2e), and we have written the integral in terms of a set of ccj 
kinematic variables that are appropriate for the all-massless pentagon; the ctj are defined by 

s M+ i = — , i = 1,2, ... ,5 (mod 5), (48) 

where s^j+i are the five independent momentum invariants for the pentagon. The coefficients 7$ 
are then 

7i = ai-2 - oti-i + on - a i+1 + a i+2 , (49) 

and the normalization constant is N 5 = 1. Plugging equation (47) into the formula (44) for the 
massless scalar pentagon, and using the dilogarithm identity Li 2 (l — x) + Li 2 (l — x~ l ) = — \ In 2 x 
to simplify the expression, we find 



/ 5 massless = r r p a]^ [1 + 2 Li 2 (l - ^i) + 2 Li 2 (l a *" 1 ^ ^ 



OLi / 6 



+ 0(e). (50) 



This solution can also be obtained by solving the differential equations (42) or (43). In terms of 
momentum invariants, the unreduced massless scalar pentagon integral is 

rmasslessrn _ r r {-Sl2Y{-S^lY 



jmassless j^j 



(-S 23 ) 1+e (-5 3 4) 1+£ (-545) 1+£ 



I + 2L l2 (l-£ 2 l)+2Li 2 (l-^)-^ 
+ cyclic permutations + 0{e). 

(51) 

jmassiess^j - g man if es tiy re al in the region where all s^- < 0; its value in physical regions can be 
obtained by the usual prescription — > s^- + ze. 

There are a few subtleties to obtaining tensor integrals via the derivative formula (30) when 
n > 4. Some of the subtleties are associated with the 1/e pole in T(n — 3 — m + 2e) for m > n — 3, 
where m is the degree of the Feynman parameter monomial being integrated. It might appear that 
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the calculation of such integrals to 0(1) would require knowledge of the scalar integral to 0(e); 
however, this is not the case. These issues are treated in more detail in ref. [2]; here we simply 
note that equation (20) for n = 4 allows the determination of h[ai\ to 0(1), given and 
to 0(1). The integrals of monomials with degree greater than 1 can then be obtained to 0(1) by 
differentiating the integrals h[ai\. Similarly, for n = 5 it is possible to derive an equation that 
yields I^aiaj] to 0(1), given 7 3 [1] and (or equivalently /f =6 [l]) to 0(1). In this case, the 
D = 6 scalar pentagon I® =6 [1] appears on the right-hand-side of the equation as well, but it is 
possible to show that for Feynman parameter polynomials that arise from tensor integrals in the 
loop- momentum, the coefficient of I® =& always vanishes. The integrals of monomials with degree 
greater than 2 can be obtained by differentiating isfojOj]. One also uses equation (42) for n = 5 
and e — > e — 1 to eliminate the derivatives dl§ =6 ~ 2e /dcti in favor of box integrals plus the scalar 
integral I® =6 ~ 2e . After doing this, the coefficient of the D = 6 scalar pentagon integral will always 
vanish, for Feynman parameter polynomials arising from tensor integrals in the loop- momentum. 

For tensor integrals with n > 6 there is another subtlety, associated with the appearance of 
A n in the denominator of the coefficient Cjj in equation (40), since A n = for D = 4 kinematics 
and n > 6. On the other hand, the scalar integrals are manifestly free of singularities as A n — > 
(so long as s^i+i 0), and by considering the tensor integrals in terms of loop-momenta, one can 
see that they also can have no singularity as A n — > 0. Therefore the 1/A n factors that appear 
in some representations of the tensor integrals (such as equation (42)), should cancel out when all 
quantities are evaluated explicitly. We have checked that this is true for insertions of up to two 
Feynman parameters. 

In conclusion, we have derived simple equations relating 7 n [l], the one-loop n-point scalar 
integral in D = 4 — 2e, to /^^[l] and the (6 — 2e)-dimensional integral I^ =6 ~ 2e [l]. In the context 
of dimensional regulation these equations may be used to recursively determine the scalar integrals 
for n > 4 as a linear combination of box integrals, up to 0(e) corrections. We also presented an 
approach to computing Feynman-parametrized tensor integrals via the differentiation of the scalar 
integral with respect to suitable combinations of the kinematic variables. This approach also leads 
to simple differential equations for scalar integrals, particularly box integrals. In reference [2] these 
general results are applied to the specific computation of box integrals with massless internal lines, 
but an arbitrary number of external masses, and to the computation of pentagon integrals with 
all lines massless and with one external mass. The latter integrals are of use in the calculation of 
one-loop contributions to amplitudes such as gg — > ggg and Z — > qqgg. 

We thank R. K. Ellis for discussions, especially regarding the cancellation of the six-dimensional 
pentagon integral from tensor integrals; and J. Vermaseren for helpful comments. 
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